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In this work we afford the statistical characterization of a linear Stochastic Volatility Model 
featuring Inverse Gamma stationary distribution for the instantaneous volatility. We detail the 
derivation of the moments of the return distribution, revealing the role of the Inverse Gamma law in 
the emergence of fat tails, and of the relevant correlation functions. We also propose a systematic 
methodology for estimating the parameters, and we describe the empirical analysis of the Standard & 
Poor 500 index daily returns, confirming the ability of the model to capture many of the established 
stylized fact as well as the scaling properties of empirical distributions over different time horizons. 
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I. INTRODUCTION 

A large number of empirical studies has shown that 
financial time series exhibit statistical features strongly 
departing from the Gaussian behavior. This finding dates 
back to the work of Mandelbrot [T] , whose attention was 
mainly focused in recognizing realizations of stable pro- 
cesses, and to the analysis of Fama [2] concerning the 
long tailed nature of the Dow Jones Industrial Average 
single components. Since these fundamental contribu- 
tions the modeling of financial returns has considerably 
grown, and very heterogeneous models, able to repro- 
duce the degree of asymmetry and the excess of kurtosis 
of the measured distributions, have been proposed. A 
non exhaustive list includes approaches developing from 
specific distributional assumptions, as it is the case of 
the Levy flights [3H5] , the Generalized Student-i or Tsal- 
lis distributions [BHH] and the exponential one [9]. Past 
empirical analysis have also proved the existence of non 
trivial scalings of higher order correlations between re- 
turns at different times, pointing toward the existence of 
a secondary stochastic process, as fundamental as that 
of the price, governing the volatility of returns. Many 
effective mechanisms allowing to reproduce the observed 
correlation structures, where the stochastic nature of the 
volatility plays a central role, were proposed. Discrete 
time models include ARCH-GARCH processes [TOl ITT] , 
and multifractal models [12j [13] , inspired by the cascades 
originally introduced by Kolmogorov in the context of 
turbulent flows. As far as continuous time approaches 
are concerned, fractional Brownian motion and stochas- 
tic volatility models have been extensively analyzed. For 
a review of the latter approach we suggest [H] and re- 
mind the reader to the discussion in Section [TTJ Focusing 
on the continuous time stochastic volatility framework, 



in this work we aim at reproducing many of the above 
mentioned facts which are generally accepted as univer- 
sal evidences, shared among different markets in different 
times. 

The structure of the paper is the following. After in- 
troducing a general class of stochastic models driving the 
evolution of the volatility, in Section |ll] we concentrate 
on a linear one able to reproduce an Inverse Gamma 
distribution in the long run. In Section III we detail 



the derivation of the moments of the probability density 
function p(x; t) of the returns over the time lag t, taking 
into account explicitly the time at which the Y process 
has started and deriving rigorously the stationary limit of 
the volatility. We describe the mechanism through which 
the power law distribution of a induces fat tails on p(x; t) 
for all the finite time lags. In Sections |IV| and |V| we de- 
rive the analytical expressions of the leverage correlation 
and the volatility autocorrelation functions respectively. 
In Section [VI] we propose a systematic methodology for 
estimating the model parameters, and we apply it to the 
time series of the daily returns of the Standard & Poor 
500 index. The relevant conclusions, along with possible 
perspectives, will be summarized in Section [VII] 



II. THE MODEL 

We consider a model where the asset price 

St = S cxp (pt + X t ) 

is a function of the stochastic centred log- return X t and fi 
is a constant drift coefficient. We assume that X t can be 
modeled with the following stochastic differential equa- 
tion (SDE) 



dX t = a t dW u 



(1) 
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where at is the instantaneous volatility of the price and 
dW\,t is the increment of a standard Wiener process. 
Since Xq = 0, from the above assumption we have that 
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(Xt) — and (\n(St/So)) = \it for all t. In the context 
of stochastic volatility models (SVMs) the instantaneous 
volatility is assumed to be a function of an underlying 
driving process Y t , i.e. a t — a(Y t ). Typically, the dy- 
namics chosen for Y t corresponds to a particular case of 
the following general multiplicative diffusion process 



dY t = (aY t + b)dt+ JcY t 2 + dY t 



:dW- 



2,t ■ 



(2) 



with suitable constraints on the parameters, in order to 
ensure the well definiteness of the process. Moreover, the 
two standard Wiener processes Wi 2 are possibly corre- 
lated 



(dW lttl dW 2tt2 ) = P 5(t 1 -t 2 )dt : 



(3) 



with p 6 [—1, 1], which is necessary to account for skew- 
ness effects and for the return-volatility correlation. For 
instance, in the Stein-Stein model [T5] [16] the volatil- 
ity is linear, a t oc Y t , and Y t follows a mean reverting 
Ornstein-Uhlenbeck dynamics corresponding to a < 0, 
b > 0, c = d = 0. Under the same Y dynamics but 
with a t oc exp (Y t ) we obtain the exponential Ornstein- 
Uhlenbeck model Q21 US]. In the Heston model [13] |2"U] 
at = \jY t and Y t evolves according to a Cox-Ingersoll- 
Ross dynamics, stemming from |2]) by taking a < 0, 
b > with c = e = 0. Finally, in the Hull- White model 
the volatility has the same functional dependence as in 
Heston, but Y t has a Log-Normal (non mean reverting) 
dynamics corresponding to a > and b = d = e = 0. 

In the Econophysics literature several studies have 
been devoted to asses the statistical properties of the 
volatility (see for instance Chapter 7 in [S] and [5T]), es- 
pecially its distribution, and it has been recognized that 
the instantaneous volatiliy, measured by suitable proxies, 
distributes in good agreement with a Log-Normal or an 
Inverse Gamma law, the best fit being obtained with the 
latter [6] which is able to better capture the heavy tail of 
the empirical distribution. None of the previously cited 
models feature an Inverse Gamma probability density 
function (PDF) for at, even though this distribution has 
been considered previously in different contexts. For in- 
stance, the Inverse Gamma was introduced in the context 
of an ARCH- like evolution of the variance in [32] , and in 
the statistical modeling of financial data the marginaliza- 
tion of Normally distributed returns conditionally on In- 
verse Gamma variance was widely exploited since it gen- 
erates generalized Students distributions (see [2"51 l2"3] L 
However, as clarified by the empirical analysis performed 
in [21] . where intra-day returns are used to estimate a 
proxy for the daily volatility, an Inverse Gamma PDF 
for of leads to an overweighting of the tail region. 

Here we afford the statistical characterization of the 
simplest linear SVM able to account for this stylized fact 
about the volatility. The process ^ has been extensively 
studied and characterized in [25] where exact solutions 
for the moments of the associated PDF have been ob- 
tained allowing to study its relaxation modes toward a 
stationary distribution, if any. In particular, when a < 



and d = e = 0, with c > 0, the process ^ has indeed an 
Inverse Gamma stationary distribution, whose support 
is [0, +00) as long as b > 0. Thereby we consider the 
following SVM 



dX t = yTcYt dWij , X a = 

dY t = (aY t + b)dt + V~cY t dW 2 , t , Y tQ = y to 



(4) 



where to < 0, yt may be a fixed constant or randomly 
sampled, and the constant factor in the expression of the 
instantaneous volatility a t — \fcY t has been added for 
later convenience. As explained in [2 5) the stationary 
PDF of a t is 



n at (tr) 



Y exp(-A/CT) 



T{v) a^ 1 



(5) 



where the shape parameter v and the scale parameter A 
are given by 

, 2a , x 26 
v = 1 and A = — — . (6) 

c V c 



III. EMERGENCE OF FAT TAILS 

A major point to be discussed before presenting a de- 
tailed derivation of our results is the different role played 
by the initial time conditions for the X and Y processes. 
Since X t represents the detrended logarithmic increment 
of the price over the time lag t, it can be directly mea- 
sured from real time series, and in a natural way we can 
assume as starting point for this process the spot time 
t = 0. On the other hand, the secondary process can 
not be observed directly but some of its statistical prop- 
erties have been measured by means of suitable proxies. 
In particular, for intra-day frequencies there is no clear 
evidence of mean reversion, that is the high frequency 
volatility is very close to its asymptotic value [55] [57] . In 
order to capture this evidence, we assume that the pro- 
cess Y, driving the returns from to t, started in the 
past at to < and we will perform the limit to — > —00 
at the end. The assumption of stationarity for the at 
process in allows also to consider the returns dX t as 
identically distributed and uncorrelated, even though not 
independent variables, by virtue of the i.i.d. property of 
the Wiener increments. 

The structure of the model Q allows to compute the 
moments of the PDF of X t at all times t recursively. 
Application of the ltd Lemma to the function X™ readily 
provides 

(X?)= l -n(n-l)cj\xr 2 Y*) ds , 

and the same Lemma proves that the correlation func- 
tions between X and Y satisfy the following differential 
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equation 



d 
dt 



(xf y«> = f 9 <x?Y t «> 



cppq (xr'rr 1 ) + \p{ P - i) C (xr 2 y/+ 2 



(7) 



where we denned Fk = ka + k(k — l)c/2, = A;6 for ev- 
ery fceN, and p,g€ N. The previous equation is a linear 
ordinary differential equation (ODE) for every p and q, 
which can be solved recursively starting from the lowest 
order of p and q [28] , and whose solution involves integra- 
tion of the moments (Y t n ) = p, n (t;t ) of the Y process. 
For every n and every time t the latter can be expressed 
as a linear superposition of exponential functions 



Hn(t; to) 



£tfj n >exp -to)]- 



(8) 



it turns out that Kj 



The explicit expressions of the coefficients in the above 
expansion can be computed as explained in [25) . and 
involves the values /Ufe(toi*o) f° r 
k = while Kq 11 ^ does not. This implies that 

whenever the constants Fj are all negative, the only term 

surviving in the limit <o — ^ — o° is ano - ^h e pro- 

cess looses every information about the distribution of 
yt a . It is worth noticing that, even though the moments 
p n {t\to) are homogeneous functions of time, when to is 
finite this is not true for the solution of Eq. ^ which 
is obtained by integration from to t, with boundary 
condition (Xq Y 9 ) = for every p > |29j . 

From the analysis of Eq. ^ it can be verified that the 
moments of X can be expressed always as a superposi- 
tion of exponential functions of the starting time of the 
volatility as follows 



= £ffj">(t)exp(- 



Fjto) . 



(9) 
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Figure 1. (Color online) Scaling as a function of to of the sec- 
ond and third moment of X at t = 1 day, for a = —16.06 yr, 
b = 0.86 yr, c = 17.84 yr and p = -0.51, \a\ /c = 0.6. Yearly 
units (1 yr = 250 trading days). 



negative (j ^ 0), (X") is finite in the stationary limit 
to — > — co, otherwise it diverges 30J indicating the emer- 
gence of fat tails in the PDF of X t . The latter case applies 
when n > v = 1 — 2 a/c, as can be checked directly from 
the definition of F n . Since F n+ i > F n , when F n > 0, 
the divergence of {X") implies the divergence of all the 
higher order moments [31] , The same condition is re- 
sponsible for the divergence of the moments p n (t) of the 
volatility for n > v (see Eq. pi) in agreement with the 
fact that the stationary distribution of the volatility ^ 
is an Inverse Gamma distribution with tail index v. Here 
we see at work a mechanism in which the power law tail 
of the stationary distribution of the volatility induces fat 
tails in the return distribution for every time lag t, and 
its scaling for large |x| is compatible with a power law 
assumption 



p(x) 



>±oo 



1 

1 1+0 



The coefficients depend on the time lag t and, more 
precisely, by virtue of the linearity of the ODEs |7]), 
they correspond to a combination of exponential terms 
weighted by polinomial functions of t. In Appendix \K\ we 
report the explicit expressions of the coefficients Hy" (t) 
for the cases n = 2 and n — 3, from which it can be 
readily verified that the skewness of the PDF converges 
to zero asymptotically for t — > +oo. A messy calculation 
would show that an analogous behavior holds for kurto- 
sis. Thus the scaling of the lowest order moments is in 
full agreement with the one of the empirical distributions 
over long time horizons |S1 |S] . When t is finite the coef- 
ficients are finite quantities themselves, and all the 
relevant information about the behavior of (X™) in the 
stationary limit of Y is retained by the io~ ex P orLeri tials 
m Eq. ||J). Two cases are possible here: if all the Fj are 



This is in agreement with empirical studies about the 
distribution of returns over daily or intra-day time 
scales [31 [SI [3"2H3"4"] , and from the previous considerations 
we are able to constraint the tail index in the following 
range 



n* < B < n* + 1 



(10) 



where n* > is the largest integer satisfying n* < v. 
As an example, in Fig. Ill it is shown the scaling of (Xf*) 
and of the absolute value of (Xf) as a function of the 
starting time of the volatility, for t — 1 day and for a 
choice of the parameters corresponding to \a\ /c = 0.6. 
For this value of the ratio the tail index of the return 
distribution is 2 < B < 3 and consequently the third 
moment of the stationary distribution of the volatility 
diverges as t becomes more and more negative, while 
(Xf) approaches its finite stationary value. 
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IV. LEVERAGE EFFECT 



correlation reads 



For the linear model Q the leverage, measuring the 
correlation between returns and volatility, can be com- 
puted exactly. Since the squared increment dX 2 pro- 
vides an estimation of the instantaneous volatility, it can 
be defined through the following function 



C{r-t) 



(dX t dX 2 +T ) 
(dX?f 



(11) 



Empirically, for arbitrary t, C{r;t) is found to be expo- 
nentially decaying for positive r and approximately zero 
otherwise, meaning that a correlation exists between past 
returns and the volatility in the future and not vice versa. 
Empirical analysis shows that it is a short range corre- 
lation; more precisely, the decay time of C{r;t) is found 
to be of approximately 69 days for U.S. stocks and even 
smaller, about 10 days, for indexes [5J. 



The numerator (111 can be rewritten as 



(dX t dX 2 +T ) = c 3 / 2 <Ci, t Y t Y 2 +T ) dt 2 , 



expressing the Wiener increment as £t dt, where Q is 
a Gaussian noise with zero mean and 1/dt variance. 
Novikov theorem [35 , 36J allows to compute the expecta- 
tion involving giving us 



{dX t dX 2 +T ) = 

2pc 2 H(r) exp (ar) (Y 2 Y t+T exp [V~cA t W 2 ( T )]) 



where we defined A t W(r) = j t +T dW s , we took into 
account the correlation structure ^ and we used the 
following expression of the functional derivative of Y 



SY, 



t + T 



5Y t+T _ 

pVcH{t) exp {ar)Y t exp [VcAW2,t(r)] 



with the Heaviside step function H(t) defined as zero if 
t < and one otherwise. The expectation f(r;t,Y) = 
(Y 2 Yt+T exp [y^AjW 2 (t)]) satisfies an integral Voltcrra 
equation of the second kind, whose derivation is detailed 
in Appendix [El and the final expression of the leverage 



2 P H(t) 



exp 



2a 



b 



a + c 



fi 2 (t)exp (a 



(12) 



which inherits the explicit dependence on t from the mo- 
ments of Y . In order to compare the previous expression 
with real data, following the discussion at the beginning 
of Section |III| we take the limit to — > — oo, so that we 
can replace ^(t) and ps(t) with their asymptotic values, 
whose general expression, valid for n < is 



^n,st 



k=X 



Ak 
F k 



(13) 



ishes, and the leverage correlation reduces to 



Substitution in Eq. ( 12 ) reveals that the first term van- 

(14) 



C(t) = - P H(t 



a(2a 



b(a + c) 

where the leverage decay time reads 

2 

2\a\-c ' 



exp 



T C = 



So, the model correctly forecasts the exponential decay 
of C(t) and its vanishing for negative correlation times. 



V. VOLATILITY AUTOCORRELATION 

The volatility autocorrelation provides an estimate of 
how much the volatility at time t + r depends on the 
value it had at time t and it is usually defined as 



A(T;t) 



(dX 2 dX 2 +T )-{dX 2 ){dX 2 +T ) 
V&x[dX 2 ] Var[dA" t 2 +T ] 



(15) 



It is a well known stylized fact [T5J [371 13H] that A decays 
with multiple time scales and in particular, it shows a 
long range memory effect, vanishing over a time scale of 
the order of few years for stock indexes. 

For the model under investigation, the volatility auto- 
correlation can be computed exactly too. Recalling again 
the Novikov theorem and the fact that 8 dW\ t/<5Ci,t = U 
the correlation entering the numerator of ( 15 ) becomes 



{dX 2 dX 2 t+T )=c 2 {Y 2 Y 2 +T ) dt 2 
+ 2 pc 5 ' 2 H{T) (Y 2 Y t+T exp [y/c A t W 2 (r)] dW ht ) dt 2 , 

but, due to the presence of dWij, the second term re- 
sults to be of order 0(dt 3 ) and therefore it can be dis- 
carded. The exact expression of the autocorrelation 
function (Y 2 Y 2 +t \ can be obtained as explained in Ap- 
pendix [Cl leaving us with 
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A(r;t) = 



exp (or) J 26 
3^4 (t) - ^2(i) 2 1 a + < 



[Ml(*)Ma(*)-M3(*)] + ex P I( a + C ) T 1 



M4(*)- 



26 



26 



where the denominator of Eq. ( |15[ ) has been approxi- 
mated with Var[dX 2 ] = c 2 [3^4(tJ^- ^(i) 2 ] m view 
of the stationary limit for Y. After replacing the mo- 
ments fj, n (t) with their asymptotic expressions (13 1 we 
end with 



A(r) = 



N ie - T 



N' 



2 a 



r/r 



(16) 



where the coefficients read 

(4a 2 - lac - 3c 



D 

N 2 =a, 



(o + c) 



(2a + 3c) (2a + c) 



and we also defined the two volatility autocorrelation 
time scales as 



1 



and 



1 



2 b 



At this point it is crucial to notice that in deriving 
Eq. (14) and Eq. (161 we assumed implicitly that the 
moments of Y t up to the order n = 4 do converge asymp- 
totically. Recalling the expression of the shape parameter 
v in (Jflj), we see this assumption imposes 



|0| 

c 



> 



(17) 



which has to be interpreted as a consistency relation for 
the model. This constraint imposes the following strict 
ordering between the time scales of the model 



< TT < -T 



with 



(18) 



where the second inequality for r-f 4 follows from the con- 
vergence of third moment of Y t which requires \a\ jc > 1. 

The expression obtained for A fails to capture the per- 
sistence of this correlation identified in several analysis 
reviewed in |39j . The lacking of power law scaling would 
not be, in principle, a serious drawback as far as one of 



the two time scales involved in ( 16 1 was sufficiently long 



However, the ordering (18 1, which is peculiar of the con- 
sidered model, makes these scales too close each other 
and the volatility autocorrelation to decay as fast as £, 
an undesired feature shared with other models, such as 
the Stein-Stein one. The persistence of A can be ac- 
counted for introducing a non linear volatility, as it is for 
the exponential Ornstein-Uhlenbeck model [TH], or cou- 
pling a third stochastic equation driving the dynamics of 



I 

the long run value of Y t as in [37] . A further possibility 
to induce a non exponential time decay would be to con- 
sider a non linear drift term for the dynamics of Y t , even 
though the analytical tractability of the present model 
will not be preserved. 



VI. ESTIMATION OF PARAMETERS 

Now we provide a systematic methodology for estimat- 
ing the model parameters, which arc the constants a, 6, 
c entering the dynamics of Y t , plus the correlation coef- 
ficient p. We perform the estimation over the Standard 
& Poor 500 (S&P500) index daily returns from 1970 to 
2010, approximating dX t with AX t — X t +At — X t 



dX t « AX t = In 



S-. 



t+At 



- ( In 



S t 



+At 



St 



where At = 1/250 yr (one trading day). Taking into 
account that aWi jt is independent of a t and that |AWi| 
is distributed accordingly to a Folded Normal law, the 
following relations hold for the model Q 



.4: 
B 

C 



IAAI 



7T 

2A7 



<AA 2 ) _ (AX 2 ) 



|AX|) = -, 
26 2 



(aw^ 2 ; 



AX\ 



At 



(2a + c)a 



(2At) 3 



AX I 



\AWi 

2b 3 c 3 / 2 
(a + c) (2a + c) a 



The constants A and B can be measured directly from the 
data, providing us an estimation of the ratio a/c through 
the relation 



D = 



B 



2 {A 2 - B) 



The value of these quantities extracted from the series 
of the daily returns of the S&P500 index are reported 
in Table U It is crucial to observe that the value ob- 
tained for the ratio \a\/c is compatible with the con- 
straint (17), supporting the consistency of our model and 
the convergence of the volatility autocorrelation. More- 
over, the same ratio provides an estimate of v = 4.579 
(see Eq. ([6])) implying for the order of the highest 
converging moment a value n* — 4; consequently, rela- 
tion ( 10 ) indicates the following range for the tail index 
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Estimators S&P500 daily returns 



A 
B 

C 
\a\/c 



0.1457 yr' 1/2 
0.0295 yr" 1 
0.0107 yr" 3/2 
1.7895 



Table I. Estimates from return sample averages. We compute 
the value of the estimators A, B, C and D for the daily log- 
returns of the S&P500 index during the period 1970-2010, 
exploiting the means of \AX\, AX 2 and |AX| 3 . 



of p(x) 



Estimators S&P500 daily returns 





0.0864 yr 


£(0+) 


-30.9515 



Table II. Estimation of the leverage time scale and its limit 
for t — > 0, obtained from the fit of the empirical leverage 
correlation (111 for the daily log-returns of the S&P500 index, 
with the model predicted expression (1141). 



Parameter Estimate from S&P500 



-16.0608 yr" 1 
0.8627 yr" 1 
8.9749 yr" 1 

-0.5089 



4 < 8 < 5. 



The leverage correlation ( 14 ) provides a way to obtain 



the two further relations needed to fix the four free pa- 
rameters of the models. Indeed, a two parameters fit of 
the function C(t) gives estimates for the time scale r 
and for the limit r — s- + 



£(0+) 



-P 



a(2a + c) 
b(a + c) 



with the results reported in Table |TT] and Fig. [2] In 
particular, the value obtained for the leverage time scale, 
t c w 21 days, and for its amplitude £(0 + ) are consistent 
with those quoted in past analysis of different stock in- 
dexes such as the Dow Jones Industrial Average [TH |37] , 
and confirm the short range nature of this effect. 

At this point all the parameters can be recovered 
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Figure 2. (Color online) Best fit of the empirical leverage 
correlation with the model prediction (141 as a function of 
the two parameters r c and £(0 + ). 



Table III. Model parameters estimated from the daily log- 
returns of the S&P500 index during 1970-2010 through the 
relations |l9|-([22l. 



through the following relations 



D 



a = cD 

a + cC 
b = — ■ 



b(a + c) 
a (2a + c) 



£(0+). 



(19) 
(20) 
(21) 

(22) 



The final results, reported in Table III show a nega- 
tive correlation coefficient, in agreement with the known 
leftward asymmetry of daily return distributions. More- 
over, our calibration provides for the relaxation time of 
the volatility process a finite value T a = — 1/a m 15 days, 
implying that, from a practical point of view, the limit 
to — > — oo is equivalent to to <C — t ct . The fitted val- 
ues of t and £(0 + ) provides a good description of real 
data, as shown in Fig. [2} on the other hand, Fig. [3] shows 
that the theoretical volatility autocorrelation for the es- 



timated values of the parameters, Eq. ( 16 1, does not cap- 



ture the long range persistence of the empirical volatility, 
as expected from the constraints (18), while it describes 



correctly the exponential decay for small values of r. 

Finally, it is important to compare the return PDF pre- 
dicted by the model with the data sample from which the 
model parameters were estimated. Since we model the 
return dynamics for increasing t, it is even more impor- 
tant to asses to which extent the diffusion process Q is 
able to capture the scaling properties of the empirical dis- 
tribution over different time horizons. At this aim, with 
the parameters fixed from the daily S&P500 series, we 
reconstruct the theoretical PDFs simulating the process 
at different time scales, t ~ 1,3, 7, 14 days, and we com- 
pare them with the corresponding empirical distributions 
obtained aggregating the daily returns. This comparison 
is shown in Fig. [4] and Fig. [5| The daily distribution is 
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Figure 3. (Color online) Theoretical prediction for the volatil- 
ity autocorrelation function of the daily returns of the S&P500 
index 1970-2010, Eq. (T6}. 
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Figure 4. (Color online) Linear plot showing the comparison 
between the return PDFs predicted by the model (lines) and 
the data for the S&P500 index, for different time scales. 

very well reproduced by the theoretical PDF, which is 
able to fully capture the leptokurtic nature of the daily 
data. The plots also confirm that the diffusive dynam- 
ics Q, once the parameters have been fixed at the daily 
scale, follows closely the evolution of the empirical curves 
for larger t. In particular, it captures the progressive 
convergence in the central region to a distribution with 
vanishing skewness and kurtosis. 

VII. CONCLUSIONS 

In this work, we have introduced a class of SVMs where 
the volatility is driven by the general process with mul- 
tiplicative noise analyzed in detail in |25j . More specifi- 
cally, we focused on the set of parameters resulting in an 
Inverse Gamma stationary distribution for the at process. 



Figure 5. (Color online) In log-linear scale, return probabili- 
ties for the model (lines) vs S&P500 returns (points) . Curves 
have been shifted for sake of readability. 



We provided an analytical characterization of the mo- 
ments of the return distribution, revealing the role played 
by the power law behavior of the Inverse Gamma in the 
emergence of fat tails. Nevertheless, even though the 
highest order moments of X diverge for every time lag, 
the analytical expressions we obtained reveal the vanish- 
ing of both the skewness and the kurtosis, in agreement 
with the Normality of returns for long horizons. As far as 
the estimation procedure is concerned, it is worth notic- 
ing that we do not exploit directly the statistical proper- 
ties of the instantaneous volatility which is an hidden pro- 
cess, but on the contrary we infer the Inverse Gamma pa- 
rameters from well established robust stylized facts hold- 
ing at the daily scale. Indeed our model correctly predicts 
zero autocorrelation for the returns, and the short range 
exponential decay of the leverage. The persistence of the 
volatility autocorrelation over yearly horizons is not cap- 
tured, and in this perspective we would like to explore 
the possibility of coupling a third SDE in the same spirit 
of [37] . Moreover, we expect that relaxing the time homo- 
geneity of the processes, as done in [35], we may induce 
time scalings more general than the exponential one. We 
also expect the proposed dynamics to be a good candi- 
date to describe the price and volatility dynamics even 
at higher frequencies. This belief is supported by the em- 
pirical analysis discussed in the literature [5j concerning 
the statistical properties of volatility proxies for intra-day 
data. A further perspective would be to explore possible 
ways to characterize analytically the PDF associated to 
the process Q or its characteristic function. This task 
requires to solve the Fokker-Planck equation for the PDF 
or its equivalent version in the Fourier space, analogously 
to what has been done in for the Heston case. Such a 
result would also allow for an application of the model in 
the context of market risk evaluation, possibly exploiting 
efficient Fourier methodologies such as those proposed 

in mm]. 



Appendix A: Coefficients of (X t 2 ) and (Xf) 

Here we report the explicit expressions of the coeffi- 
cients Hj 71 ^ (t) entering the expansion ^ of the moments 
of X t for the cases n = 2 and n = 3. They were used to 
plot the analytical curves in Fig. [T] 

H< 2) (*) = cK ( 2 h 
H[ 2 Ht) = cK[ 2 > [ °*PW)-1 



r exp(F 2 t) - l j 
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+A 2 
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where the coefficients and Kj , entering the ex- 

pansion {[sj) of the moments of Y t , read 
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in the form 

f(T,t;Y) = (y 2 (y t + J* +T (aY s + b) ds^j cxp [VcA t W 2 (r)]\ 
+ (y 2 (s/i J* +T Y s dW2, s ~J cxp [VcA t Hf 3 (T)]\ 

= (cxp [VSA t W 2 (T)]> [ M3 (t) + 6TM2(t)] 

+ a (y 2 Y s cxp [v^A t W 2 (r)]^ ds 

(Y t 2 Y s cxp [V^AtW 2 (r)] dW 2>s ) . 

Taking into account that for t < s < t + r we can always 
split A t W2(T) as 

A t W 2 (r) = W 2 , t + T - W 2 ,« = W 2 , t + T - W 2>s + W 2iS - W 2 , t 
= A 3 W 2 {t + t - s) + A t W 2 (s - t) , 

the function /(r, t;Y) becomes 

f(r,t:Y) = (cxp [^A t W 2 (r)]) + 6r M2 (t)] 

+ <V t 2 y t+T ,cxp [v^A t W 2 (r')]) 

X (cxp [v / cA t + T /H , 2 (r - r')] ) dr' 
+ V^f {Y 2 Y t + T , cxp [v^A t lV 2 (r')]) 
x (cxp [^A t + T ,W 2 (r - r')] <JW 3)t+T /) , (Bl) 

where we changed the variable of integrations to r' = 
s — t. Since the process y/c A t+T 'W 2 {T — t') is Normally 
distributed with zero mean and variance c(t — t'), and re- 
calling the expression of the Gaussian characteristic func- 
tion, cf> , we can write 

(exp[ N /5A t+T ,VK 2 (r-T')]>= G M \^ = _. = cxp [j (r - r')] . 

Application of the Novikov theorem also gives 



/<5cxp [v^ C2, s dsl 
(cxp [VHA t+T , W 2 ( T - r')] ^ 2 , t+r , ) = ^ , L Cw2(t + ; T0 

= Vc cxp [|(t - r')j , 

where we expressed the Wiener variation in terms of a 
Gaussian white noise £2,4 as rfW~2 ,t = £ 2 .t cfa. Replacing 
the previous expressions in Eq. (Bl) we conclude that 
/(t, t; Y) has to satisfy 

/(r, t; y) - (o + c) ^ /(t', t; y) cxp [j(r - r')j efr' = 

|> 3 (i)+6rf* 2 W] , 



which is a Volterra equation of the second kind, whose 



solution leads to Eq. ( 12 ) 



Appendix C: Computation of (YtYt+ T )- 



Appendix B: Derivation of Eq. (12 1 



After expressing Y t + T in terms of its integral solution 
form t to t + t, the function /(r, t; Y) can be rewritten 



With reference to the model the cross correlation 
(F t m F t ^L r ) can be computed exactly. Provided to express 
Yf+T as integral solution from t to t + t 



A„y™ + A n Y™~ 



1 



. dW 2 . 



it is straightforward to check that (Y^ n Y^i rT ) satisfies 
the following equation 



dr 



(CI) 



whose substitution in Eq. (CI) provides the solution 



(Y?Y? + \= C * V (F 2 T)m{t) + [cxp (ar)- cxp (F 2 t)] M3 (*) 

\ > a, — 1*2 



-{tt [cxp(F 2 r) - 1] — [cxp (ar)-cxp (F 2 t)\ 



F 2 



a - F 2 



} 



which is an ODE provided that the correlation 
^my^n-i^ nag been computed at the lower order n — 1. 
In particular, for the case m = n = 2, we need the fol- 
lowing correlation 



(Y^Y t+T ^ = cxp {ar)fj, 3 (t) [1 - cxp (or)] /J 2 (t) ; 
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